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Abstract 

By means of the Nyquist method, we investigate the linear stabil- 
ity of electrostatic waves in homogeneous equilibria of quantum plas- 
mas described by the Wigner-Poisson system. We show that, unlike 
the classical Vlasov-Poisson system, the Wigner-Poisson case does not 
necessarily possess a Penrose functional determining its linear stabil- 
ity properties. The Nyquist method is then applied to a two-stream 
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distribution, for which we obtain an exact, necessary and sufficient 
condition for linear stability, as well as to a bump-in-tail equilibrium. 



PACS numbers: 52.30. -q, 52.35. -g, 52.90,+z, 05.60.Gg 



1 Introduction 

The topic of quantum plasmas has recently attracted considerable attention 
0J-H. A central reason for this accrued interest derives from the importance 
of quantum effects in the performance of today's micro-electronic devices, 
for which classical transport models are not always adequate in view of the 
increasing miniaturization level that is now entering the submicron domain. 
Hence, it is desirable to achieve a good understanding of the basic prop- 



erties of quantum transport models. The Wigner-Poisson system p~Q[]- p~2 
is a quantum transport model that has proven to be suitable in the treat- 
ment of quantum devices like the resonant tunneling diode |T|. Moreover, 



it has been referred [|13j to as perhaps the only kinetic quantum transport 
model amenable to detailed numerical simulation. In the present work, we 
address the question of the stability of small-amplitude waves, described by 
the Wigner-Poisson system. 

A convenient tool to investigate the linear stability of systems having a 



dispersion relation is provided by the Nyquist method |14|],[[U|]. Let us briefly 
review the basis of this approach. Let D(uj,k) = be the dispersion rela- 
tion, where uj and k are the frequency and wave-number for small-amplitude 
oscillations. In most practical cases, it is impossible to solve exactly the dis- 
persion relation for uj as a function of k, some kind of approximation being 
necessary. Hence, the imaginary part of the frequency, which determines the 
stability properties of the system, can be obtained only in an approximate 
way. However, exact results can be found by splitting D in its real and imag- 
inary parts, D(uj, k) = D r (uj, k) + iD { (u, k). Then, for fixed k and real uj, by 
varying uj from minus to plus infinity we can draw a diagram in the D r x Di 
plane. The resulting curve, known as the Nyquist diagram, determines the 
number of unstable modes of the system, which equals the number of times 



the origin is encircled by the diagram |14| . For example, using the Nyquist 
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method, one can show that equilibrium distributions that are monotonically 
decreasing functions of the energy are stable against small perturbations. 
Moreover, for symmetric equilibria with at most two maxima, the sign of the 
so-called Penrose functional ||14|| , |15|1 determines the linear stability of the 
classical Vlasov-Poisson system. 

In view of the utility of Nyquist's method for classical plasmas, it seems 
desirable to investigate whether it can be applied to the quantum case too. 
This approach is justified, since the linear stability of waves in the Wigner- 
Poisson system is described by a dispersion relation, and is therefore amenable 
to Nyquist's treatment. However, we cannot a priori expect to obtain a result 
as general as in the classical case. Indeed, as we shall see, the question of 
stability is subtler in the quantum framework, a typical example being pro- 
vided by the two-stream instability [[16]. For simplicity, in the present work 
we shall only consider homogeneous equilibria for one-dimensional electro- 
static plasmas consisting of mobile electrons. An immobile ionic background 
guarantees overall charge neutrality. 

The paper is organized as follows. In Section II, we develop the funda- 
mentals of the Nyquist method as applied to quantum plasmas described 
by the Wigner-Poisson system. The stability properties of quantum plas- 
mas are determined by the specific form of the quantum dispersion relation 
JT7| ] , [Tj| . We show that there are a rich variety of possible behaviors in quan- 
tum systems, which are not present in classical Vlasov-Poisson plasmas. In 
particular, in Section III, we prove that a quantum analogue of the Penrose 
functional cannot exist. To show this, we consider symmetric equilibria with 
at most two maxima. Nevertheless, the Nyquist method can still be used 
for Wigner-Poisson plasmas. This is explicitly shown in Section IV, where 
we study a two-stream equilibrium, described by a bi-Lorentzian distribu- 
tion function, which is amenable to exact calculations. We find an exact 
criterion for stability, which reduces to the classical criterion when quantum 
effects becomes negligible. However, large quantum effects can destroy the 
instability occurring in the purely classical case. In Section IV, we also in- 
clude the example of the physically relevant distribution corresponding to a 
bump-in-tail equilibrium. Our conclusions are given in Section V. 
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2 Quantum dispersion relations 



If f(x,v,t) is the Wigner quasi-distribution and the scalar potential, then 
the Wigner- Poisson system [|Kj-|]l2| reads 

' ) f + v U: = [ dv'K{v'-v,x,t)f{v',x,t), (1) 

dvf-no), (2) 



dt dx 

d 2 d> e 



dx 2 Eq 

where K(y' — v, x, t) is a functional of the scalar potential, 

im(v' — v)X 



K(v' — v, x, t) 




H J 

~,*)-<Kz + ^)J • (3) 



Here, no is a background ionic density, — e and m are the electron charge 
and mass, h is the scaled Planck constant and Eq is the vacuum dielectric 
constant. We take periodic boundary conditions in space and assume that 
for large \v |, / and all its velocity derivatives tend to zero. We also assume 
that the initial Wigner function is everywhere positive. However, the time 
evolution determined by Eq. (Jl|) may force the Wigner function to assume 
negative values. Hence, a strict interpretation of / as a true probability 
distribution is impossible. In spite of that, the Wigner function may be used 
as a useful mathematical tool to compute macroscopic quantities such as the 
charge density and electric current. 

The linear stability of a plasma, be it classical or quantum, is determined 
by the dispersion relation, which is obtained after Fourier transforming in 
space and Laplace transforming in time. Following this procedure, we obtain 



T7| for a frequency lu and a wave number k 

D(k,u) = D r (k,u)+iD i (k,u) = 0, (4) 
where the dispersion function D(k,u) is given by 

^ tt n t u l f dvFiv) 

D r (k,u) = l-^W t /M2 L 2IA 2 . 5 

n k z Jp (v — uj/ky — h k 2 /4m 2 

f i \ ( -nfU hk, ^,uj %k.\ 
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In Eq. (|5[), P stands for the principal value symbol and F(v) denotes the (spa- 
tially homogeneous) equilibrium Wigner function. Also, cu p = (noe 2 /me ) 1 ^ 2 
is the usual plasma frequency 

The quantum formulae reduce to the classical ones as h — > 0. In particu- 
lar, 

Moreover, no matter what the value of h, for \u\ — > oo we have D r — > 1 and 
— > 0, as in the classical case. 

The topology of the Nyquist diagram is determined by the sign of D r at 
the points where Di = 0. As mentioned in the Introduction, the number of 
unstable modes equals the number of times the Nyquist curve encircles the 
origin. Therefore, unstable modes can only exist if D r < for at least one of 
the points where Di = 0. In the classical case, the zeroes of the imaginary 
part of the dispersion function are determined by the points at which the 
distribution function has zero derivative. In the quantum case, according to 
Eq. (HD, the decisive points are the real roots v of 

F(v + H) = F(v - H) . (8) 

Here and in the following, 

The geometrical interpretation of Eq. @ is simple: we have to find the 
points Vq that are equidistant to any two points at which F has the same 
value (see Fig. [I]). The corresponding distance is H. In a sense, Eq. (||) 
is the finite difference version of the classical condition dF/dv(v = vq) = 0. 
Finally, as Nyquist's diagram is obtained taking exclusively real frequencies, 
only the real roots of Eq. (§) are relevant. 

The basic tasks we have to perform are first solving Eq. (|8p, obtaining 
all real roots vq for a given H, and then studying the sign of D r at each such 
root, taking u = kvo. Using Eq. (H), we have 



J*L[ dv EM 

n k 2 Jp (v — Vq) 2 — H 2 



D r (k,u = kv ) = l v - / dv-. y J — . (10) 
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Now, in the Cauchy principal value sense, 



dv 

p (v- v ) 2 -H 2= °' (H) 

Using this fact, we can rewrite the real part of the dispersion function in the 
more convenient way 



D r (k, u = kv Q ) = 1 + —75 / dv ; ° " „\ "> . (12) 



^ r dv F(v + H)-F{v ) 
n k 2 J (v — v ) 2 — H 2 

In this form, the principal value symbol is not needed anymore, since the 
integrand is regular as v goes to v ± H. Indeed, using the fact that F(vq + 
H) = F(v - H) from Eq. (|), we find that 

Fto + H)-F(«) = 1 dF 
v^v ±h ( v - Vq )2-H 2 2H dv y ' v ; 

is a finite quantity. A similar (but not identical) regularization procedure 
holds in the classical case too [fH| . 

Equations (||) and ( [X2j ) are the fundamental equations for Nyquist's method 
for one-dimensional quantum plasmas, in which only electrostatic fields are 
present. So far, the treatment has been completely general. Let us now con- 
sider some particular equilibria in order to analyze the consequences of Eqs. 
<E) and (HI. 



3 Equilibria with one or two maxima 

If the equilibrium Wigner function F(v) has a single maximum v max , then 
the geometric meaning of Vq is sufficient to show that Eq. (|8]) has always one, 
and only one, real solution v for any value of H (see Fig. |1|). Depending 
on the shape of F, this solution can differ considerably from v max (one has 
^0 — v max when F is symmetric with respect to v max ). However, as H goes 
to zero, and again from geometrical arguments, we can convince ourselves 
that Vq approaches v max . Indeed by definition u is equidistant to the points 
v' and v" for which F(v r ) = F(v"). The corresponding distance from vq to 
either v' or v" is H . 
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Furthermore, for (v — v ) 2 > H 2 we have F(vq + H) > F(y) and for 
(v — Vq) 2 < H 2 we have F(vq + H) < F(v). Hence, the integrand in Eq. 



Since there is only one root for Eq. (f|), we deduce that the Nyquist diagram 
cannot encircle the origin, and therefore no unstable modes can exist for an 
equilibrium with a single maximum. Thus, no matter how strong quantum 
effects are, the conclusion is the same as for the classical case. 

Let us now consider equilibria with a single minimum, v m i n (see Fig. 0). 
This is equivalent to consider equilibrium Wigner functions with only two 
maxima, as on physical grounds the equilibrium function must decay to zero 
as \v\ — > oo. Physically, such equilibria correspond to a situation where two 
counterstreaming electron populations (with similar temperatures) co-exist. 
In the classical case, the Nyquist diagram for this situation leads to the 
construction of the so-called Penrose functional 



which determines the stability properties of the system. The inequality 
P[F] < is a necessary and sufficient condition for instability, for appro- 
priate wave numbers. This can be easily seen from the classical limit H — > 
of Eq. (|i"2D. Classically, the points Vq where Di = are the maxima (±v max ) 
and the minimum (f m j n ) of the equilibrium distribution. For vq = ±v max , 
the integrand in Eq. ( |12|) is always positive, and thus cannot lead to in- 
stability. For vq = v min , the real part of the dispersion function reduces to 
D r = 1 + (u 2 /nok 2 ) P[F}. If the Penrose functional is positive, instability 
is ruled out. If it is negative, one can always choose k small enough so that 
D r < and therefore some unstable modes must exist. This completes the 
proof of the necessary and sufficient Penrose criterion. 

The natural question now is whether there exists an analogue Penrose 
functional for the quantum case. For simplicity, in the following we restrict 
our discussion to Wigner equilibria that are symmetric about v min , the point 
at which F attains its minimum value. By a Galilean transformation, this 
point can be taken as v m i n = without loss of generality. We first notice that, 
in the classical case, one only has to consider the three velocities for which 



flT2| ) is always positive, implying that the real part of the dispersion function 
is a positive quantity. Also, for \u\ —>■ oo we have D r —>■ 1 and Di — > 0. 




(14) 
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the equilibrium distribution function has zero derivative. In the quantum 
case, however, depending on the shape of the equilibrium Wigner function, 
there can be more then three roots for Eq. (Bp, with fixed H. For instance, 
in Fig. |2|, root v\ (connecting one increasing and one decreasing branch of 
the distribution) can be obtained from the local maximum v max , by varying 
H continuously from zero to a certain value. The root t>2 (connecting two 
decreasing branches of the distribution) is of a different nature, arising only 
for sufficiently large H. Indeed, it is not difficult to realize that, in the case 
of two maxima, there are always only three roots for Eq. (|8]) if H is small 
enough, and up to five roots for larger values of H . Also notice that, for 
symmetric equilibria, the point v = is always a root, irrespective of the 
value of H. For a given H, possessing one, three or five real roots depends 
on the details of the equilibrium. It is not difficult to prove that, in the case 
of a two-humped distribution, a sufficient (but not necessary) condition for 
having five roots to Eq. (§) is that F{y min ) = 0. This can be shown by 
plotting the left- and right-hand sides of Eq. (H) as a function of v , and 
looking at the intersections of the two curves. In general, we obtain that 
one can have five solutions when F(v m i n ) is smaller than a certain threshold. 
Note however that five roots only appear for sufficiently large values of H; for 
small enough H, there are always only three roots. In Section IV, we shall 
examine a bi-Lorentzian distribution possessing at most three solutions. In 
addition, we shall discuss another two equilibria, which possess five solutions 
for sufficiently large H. 

Let us now consider the question of the existence of a quantum Penrose 
functional. We need to examine the sign of D r at the different solutions of 
Eq. (PJ). The root t> = v min = always exists and can yield either a positive 
or a negative value for the integral in Eq. (0), depending on the shape of the 
equilibrium and the value of H. One can actually prove that the integral can 
be negative only if H < v*, where v* is the positive solution of the equation 
F(0) = F(v*). 

We now analyze the other roots of Eq. (|8|). Let v± be the root obtained 
from the maximum of F at the right of v = by varying continuously H 
from zero to some particular value (see Fig. |2|). Referring to Fig. ^| and to 
Eq. fll2|) (with Vq — Vi), we conclude that the integrand in D r is negative for 
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—v\ — H < v < —Vi + H. Thus, in principle, the real part of the dispersion 
function can be negative. However, one could imagine that the negative 
contribution for —V\ — H < v < —V\ + H is compensated by a positive 
contribution corresponding to V\ — H < v < v\ + H. Let us examine this 
possibility. Using the fact that F is even, we obtain 



dv Fjv-y + H)- F(v) r^+ H dv F(vi + H) - F(v) 



V!-H (v — Vi) 2 — H 2 Jvi-H (v — Vi) 2 — H 2 



(15) 



vi—H (v 2 -(H + v x ) 2 ) (v 2 -{H- Vl ) 2 ) 



For vi — H < v < vi + H, we have F(v) > F{v\ + H), v 2 > (H — v{) 2 and 
v 2 < (H + Vi) 2 . Hence, the integrand in Eq. (|15|), which can give the only 
negative contribution for D r , is negative provided 

v 2 < H 2 - v\ , (16) 

which is impossible in the prescribed range of velocities, as Vi > H by con- 
struction. Therefore, we always have D r (k,u = kvi) > 0, where Vi is the 
(semi-classical) root for Eq. (§) obtained from the positive maximum of F, 
and the same argument holds for the symmetric root —Vi. This is analo- 
gous to the classical result shown above, according to which D r is positive at 
the two maxima of F(v). Indeed, the roots ±t>i coincide with ±v max when 
H -> 0. 

However, this is not the end of the story for the quantum case. Indeed, for 
sufficiently large values of H, it is possible to access the roots ±v 2 (connecting 
two decreasing branches of the distribution) shown on Fig. |^. [This is not 
in contradiction with the above statement that some equilibria only display 
three solutions to Eq. (]8|). Solutions of the type t> 2 always exist, although 
they may correspond to different values of H than vi, so that for a fixed H 
there are indeed only three roots]. For the roots ±v 2 , which are of a strictly 
quantum nature, we cannot anymore obtain, a priori, D r > 0. For instance, 
for the particular choice of v? shown on Fig. |2|, the region — v 2 — H < v < 
—V2 + H contributes a negative value to D r (k, u = kv 2 )- The same is true for 
the root —v 2 - This is because, over most of the region —v 2 — H < v < —v 2 +H 
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one has F(v 2 + H) < F(v) and (v — V2) 2 > H 2 , implying that the integral 
in Eq. flT2|) is negative. Another choice of v 2 may have yielded the opposite 
result, so that the sign of D r (k,u = kv 2 ) cannot be determined a priori. 
As the parameter H depends on the wave number, it is always possible 
to choose k so as to access a root of the type ±t> 2 , for which the sign of 
D r is undetermined. The conclusion is that there is no quantum Penrose 
functional, since the topology of the Nyquist diagram can be changed, in an 
essential way, by the value of D r at the quantum roots for Eq. @. Each 
specific equilibrium must be studied in detail. In the following Section, we 
shall illustrate the previous theory using some concrete examples. 



4 Examples of two-stream and bump-in-tail 
equilibria 

Let us consider a two-humped equilibrium given by 

F ^ = {(v-a) 2 + A 2 + (v + a) 2 + A 2 ) ' ^ 

where A is a measure of the dispersion of the distribution and a is a parameter 
associated to the distance between the two possible maxima. If a 2 < A 2 /3 
this bi-Lorentzian distribution degenerates into a one-humped equilibrium, 
which is consequently stable against linear perturbations, both in classical 
and quantum cases. The major advantage of dealing with Eq. QT7| ) is that it is 
amenable to exact calculations, thus providing an appropriate example of the 
use of the Nyquist method for quantum plasmas. Moreover, it models the 
physically relevant situation of two counterstreaming electron populations 
that co-exist within the same plasma. 

Inserting Eq. ([17]) into Eq. (|8[), we obtain the following solutions, 

v° = 0, (18) 



v, 



I = ± (H 2 - a 2 - A 2 + 2Va 2 + AVa 2 - H 2 ) 1/2 , (19) 



± (H 2 -a 2 - A 2 - 2Va 2 + AVa 2 - H 2 ) 1/2 . (20) 
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It is easy to check that the roots given in Eq. ([21]) are always complex, 
whatever the values of H, a and A. However, the roots ([TJ^) can be real, 
provided 

a 2 > \ H 2 _ A 2) + 2 (A 4 + ^2 A 2 + ^4)1/2 ^ 

3 3 



Thus, there can be one or three relevant roots, according to condition (|2lD. 
This inequality, when satisfied, can also be seen to imply a 2 > A 2 /3, which is 
the same as the condition for the existence of two maxima. Hence, there can 
be three real roots if and only if F is two-humped, which is not surprising in 
view of the arguments given in the preceding Section. 

An equivalent and illuminating way to rewrite Eq. (|21|) is 



H z < v z max , (22) 
where v max denotes the (positive) point where F is maximum, 



a 



+ A 2 ) 1 / 4 (2a - Va 2 + A 2 ) . (23) 



1/2 



Hence, 2H cannot exceed the distance between the two maxima of F. Notice 
that the right hand side of Eq. (^) can be real only if a 2 > A 2 /3, that 
is, if there are two maxima, which is again a natural result. For very large 
quantum effects, only the root v ° = survives. 

As there is no quantum Penrose functional, it is necessary to calculate 
D r at all possible roots (|l8l)-(|l^). We obtain 

DAk, „ = kv° = 0) = 1 + g w _ Jf + * { ^% + A , ■ (24) 
which can be negative if and only if 

a 2 > A 2 + H 2 . (25) 

In addition, 



D r (k, u = ±kvl) = 1 + 16 ^ '— , (26) 
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where 5 and U 8 are the positive-definite quantities 

S = Va 2 + A 2 - Va 2 - H 2 , (27) 
U* = (((v -H) 2 -a 2 ) 2 + 2A 2 ((v -H) 2 + a 2 )+A 4 )x 

x ({{v + H) 2 - a 2 ) 2 + 2A 2 ((^ + H) 2 + a 2 ) + A 4 ) . (28) 



We can show that Eq. fl2ll) implies a 2 > H 2 , so that D r as given by Eq. (|26|) 
is indeed always positive. 

In view of Eqs. (p4j) and (0), we see that Eq. (|25| ) is a necessary and 
sufficient condition for linear instability. That this condition is sufficient can 
be easily proven: suppose that we have found a wave number k satisfying 
Eq. (^5|); then any k < k will also satisfy it; by taking k small enough, we 
can make the second addendum of Eq. fl24|) (which is negative) arbitrarily 
large in absolute value and therefore obtain D r < 0. Note, however, that 
putting an equality sign in Eq. (25) and solving for k does not provide the 
transition wave number between stable and unstable behavior. In order to 
obtain it, one has to set Eq. (|24]) to zero and solve for k. 

Equation ( p5|) means that the plasma can become unstable for sufficiently 
large a (the two maxima are sufficiently far apart), small A (small disper- 
sion) or small H (small quantum effects). We also notice that, as H depends 
on the wave number, quantum effects can suppress the instability for small 
wavelengths. The instability condition Eq. (|25"D confirms the numerical re- 
sults by Suh et al. . Here, however, we have derived an exact analytical 
criterion for quantum linear stability of a two-stream equilibrium. 

On Fig. H|, we have plotted the Nyquist diagrams for the two-stream 
equilibrium of Eq. (|T7|) with a = 3, A = 1, k = 0.2 and four different values 
of h (units for which e = Eo = m = n = l are used). We observe that 
stabilization of the k = 0.2 mode occurs somewhere between % = 25 and 
h = 27. This is in agreement with the previous formulae: indeed, with this 
set of parameters, it is found that D r (k,u = kv$) changes sign for h ~ 25.5. 
Furthermore, Fig. [3] also shows a change in the topology of the Nyquist dia- 
gram. For figures (a) — (c), the diagram intersects the horizontal axis in three 
points [excluding the point (1,0) that corresponds to to = ±oo]; note that two 
such points coincide, because of the symmetry of the distribution. For figure 
(d), only one intersection survives. This change in topology corresponds to 
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having one or three solutions to Eq. (^j), which is determined by Eq. (PT|). 
The transition is found to occur for ft ~ 30, which is in agreement with the 
diagrams of Fig. ^. 

Finally, we point out that large quantum effects are not necessarily sta- 
bilizing. For the two-stream equilibrium of Eq. fll7|), with a = 3 and A = 1, 
the wave number k = 0.287 is classically stable. However, increasing h, one 
finds that D r (k,u> = hv®) becomes negative on an interval approximately 
given by 7 < % < 11.8. This can be easily verified by plotting Eq. (|24|) as 
a function of h or by direct substitution of the above values. However, this 
destabilizing effect occurs for rather limited ranges of h and k. For example, 
wave numbers k < 0.28 are classically unstable and are stabilized for large 
enough h (as in our previous example with k = 0.2); on the other hand, wave 
numbers k > 0.29 are classically stable and remain stable for any value of 
h. Only wave numbers very close to the value k = 0.287 display the unusual 
behavior described above. For this reason, we can still conclude that the 
most likely outcome of quantum effects is stabilization. 

We now show that we can explicitly write a distribution function for 
which there can exist five real roots for Eq. (^). Consider the two- humped 
equilibrium 

F( w ) = J^ 2 exp(-*; 2 /a 2 ), (29) 

where a is a parameter related to the equilibrium temperature. As F(0) = 0, 
we should expect that Eq. (j§) possesses five real roots for some (large enough) 
H. We now give an explicit proof of this fact for the above equilibrium. The 
solutions to Eq. (|8|) are obtained in this case from the equation 

tanh = G(v; H) , (30) 

where we have defined 



G(v; H) = 4^7T, ■ (31) 



v 



+ H 2 



Apart from the obvious root v — 0, we can have two or four additional real 
roots. By plotting the left- and right-hand side of Eq. (RDf) as a function of 



v (see Fig. ^), we can show that there will be a total of five real roots if and 
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only if 

^(, = 0)<A t anh(^W0). (32) 
av av \ a z J 

This implies that there can be five real roots provided H > a, that is, for 
sufficiently large quantum effects. Otherwise, only three (semi-classical) so- 
lutions are possible. 

In the remaining part of this Section, we address the question of the 
quantum linear stability of an equilibrium characterized by a large central 
distribution of electrons with in addition a small bump in the tail. This is 
a standard problem in plasma physics, with the small perturbation to the 
central distribution representing a beam injected in the plasma. Here, we 
consider the quantum aspects of the problem, by using Nyquist diagrams. 
The so-called bump-in-tail equilibrium has a single minimum, but, as there 
is no quantum Penrose functional, we are lead to compute the real part of 
the dispersion function at all critical points (zeroes of the imaginary part of 
the dispersion function). Nevertheless, the Nyquist technique is less expen- 
sive than, for instance, direct calculation of the dispersion relation, since it 
requires the value of D r at a few points only. 

To model the bump-in-tail equilibrium, we use the following distribution 
(see Fig. ||) 

2n [1 - y/2>/a)] 2 

3vra [l + iv/af] 2 ' K ' 

where a > is a reference velocity that can be scaled to unity without loss 
of generality. Henceforth, we set a = 1. The distribution of Eq. (|33]) is a 
particular case of a one-parameter family of bump-in-tail equilibria whose 
classical linear stability properties have been recently studied via Nyquist 
diagrams [pH[] . In the quantum case, there is no Penrose functional and the 
analysis is more involved. 

Inserting Eq. (^) into the determining equation Eq. (|8|), we obtain the 
following equation for v 



1 - V2(v + H) 1- y/2(v - H) 



(34) 



l + (v a + Hf l + (v -H) 2 ' 

The plus sign yields the second-degree equation : v 2 — V2v — H 2 = 0, with 
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solutions 

v±(H) = - (y/2 ± V6 + AH 2 ) . (35) 

Note that, in the limit H — > 0, these solutions correspond to the two maxima 
of the equilibrium distribution. Taking the minus sign in Eq. ([53) yields the 
third-degree equation 

V2v% - vl + y/2(l - H 2 )v - (1 + H 2 ) = . (36) 



It is easy to prove that Eq. ( |36"D has one real solution for H < 3 and three 
real solutions for H > 3. Furthermore, the largest of such solutions is always 
positive, and coincides with the minimum of the equilibrium distribution 
when H = : we shall call this solution v m (H). The other two solutions 
(which are real only when H > 3) have no classical counterpart, and will be 



called v q i(H) and v q2 (H). A graph of all the roots of Eq. flSHf ) as a function 



of H is provided on Fig. ^. Again, the existence of five real roots for some 
values of if is a consequence of the fact that F{v m i n ) = 0. 

As further calculations are rather cumbersome, we only report here the 
most relevant results (mostly obtained using the mathematical package MAPLE). 
For H > 3, we have obtained numerically, using Eq. flT2|), that D r (k,u) = 
kvgi) > and D r (k,u = kv q <i) > 0. At least for this particular example, this 
can be shown to imply that the purely quantum solutions are irrelevant to 
the linear stability properties of the equilibrium. Moreover, we have found 
numerically that D r [k,u = kv-(H)} > 0. After an involved analysis to de- 
termine the ordering of all the solutions for Eq. (§), the conclusion is that 
the unstable modes satisfy 

D r [k,u = kv m {H)) < 0, D r [k, to = kv + (H)) > . (37) 

This holds whatever the value of H. Remembering the dependence of H 
on the wave- number and taking into account the explicit forms of v m (H) 
and v + (H), it appears that the pair of conditions ( |3~?|) are very complicated 
expressions of k and ft. However, using appropriated units in which n = 
uj p = m = 1, we were able to solve Eq. (|3?p numerically for a few values 
of h (measured in units of ma 2 /uj p ). For Ti = 0, we found that the unstable 
modes satisfy 0.18 < k < 0.82, where k is measured in units of uj p /a. This 



15 



is the classical condition for linear instability. For % = 10, the instability 
range is given by 0.20 < k < 0.36. We see that the total band of instability 
becomes smaller for a non-zero Planck's constant. Further increasing h, 
taking h = 100, we found that the unstable linear waves must satisfy 0.15 < 
k < 0.17. For even larger quantum effects, there is virtually a suppression 
of all unstable modes. This is again in agreement with the numerical results 



of Suh et al. [19], where large quantum effects were shown to stabilize all 



classically unstable modes for a two-stream equilibrium. 



5 Conclusion 

In this paper, we have discussed the Nyquist method for the study of the 
linear stability of spatially homogeneous quantum plasmas described by the 
Wigner-Poisson system. For classical Vlasov-Poisson plasmas, this method 
provides a simple way to analyze the stability properties. Furthermore, for 
the special case of two-stream equilibria, one can construct a simple func- 
tional (known as Penrose functional), whose sign determines whether unsta- 
ble modes exist. 

The main conclusion of the present work is that the stability analysis of 
quantum plasmas is generally subtler than in the classical case. In particu- 
lar, we have shown that no simple analogue of the Penrose functional can be 
constructed in order to determine the stability properties of a two-humped 
equilibrium. Hence, a detailed analysis is necessary for each particular case, 
with generic and universal conclusions being more difficult to obtain. How- 
ever, we were able to prove that one- humped equilibria (i.e., with a distri- 
bution that is a monotonically decreasing function of the energy) are always 
stable: this is the same result as for classical plasmas. 

The main mathematical reason for the subtler behavior of quantum plas- 
mas is that the wave number now enters the real part of the dispersion 
function through the parameter H = hk/2m. This can change the topology 
of the Nyquist diagram, not only by varying h, but also by varying the wave 
number at fixed h. Physically, this means that new unstable modes can arise 
by resonant interaction between the quantum velocity H and some other 
typical plasma velocity. Indeed, such purely quantum unstable modes have 
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been observed |TB| for the special case of two-stream equilibrium 



Hv) = ^5(v -a) + ^6(v + a) , (38) 
where 5 is the Dirac delta function and ±a the velocities of each stream. 



This equilibrium can be amenable to exact calculation [jig] . 

Even when general or exact results cannot be obtained, the Nyquist tech- 
nique can be successfully used for the study of particular equilibria, as was 
shown in Section IV. The bi-Lorentzian equilibrium treated in that Section 
has shown that large quantum effects generally contribute to stabilize per- 



turbations [0. This is not always the case, however, and we have produced 
an explicit example of a wave number that is classically stable and becomes 
unstable for finite h. Moreover, the Nyquist method has enabled us to derive 
an exact stability criterion for such a bi-Lorentzian equilibrium. The Nyquist 
technique was also applied to a classically unstable bump-in-tail equilibrium. 
Again, large quantum effects were shown to reduce the range of unstable 
wave numbers. 
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Figure 1: Graphical representation of the geometric meaning of vq [solution 
of Eq. (Bp] for a one-humped distribution function. The distance between v 
and both v' and v" is equal to H. The Wigner function is represented on the 
vertical axis and the velocities on the horizontal axis. 



Figure 2: Semi-classical (vq = ±t>i, solid horizontal lines) and purely quan- 
tum (v = ±t> 2 , dashed horizontal line) solutions of Eq. (||) for a symmetrical 
two-stream equilibrium. Also note that v = (dotted line) is always a so- 
lution. Units are conveniently rescaled. 



Figure 3: Nyquist diagrams for the two-stream equilibrium of Eq. (|17D with 
a = 3, A = 1, k = 0.2 and h = 0.001 (a), 25 (b), 27 (c) and 40 (d) (units for 
which e = e = m = n = l are used). Diagrams (c) and (d) indicate that 
quantum effects have suppressed the instability. 



Figure 4: Plot of the left-hand side (solid line) and right-hand side (dashed 
line) of Eq. ( ]30"D as a function of v, for a = 1 and H = 0.7 (a) and H = 1.2 
(b). The inset is a zoom of the region 0.9 < v < 1.4 for case (b), showing in 
detail the extra solutions arising for H > a. 



Figure 5: Velocity distribution corresponding to the bump- in-tail equilib- 
rium of Eq. fl3"5| ) for a = n = 1. 



Figure 6: Plot of the roots v of Eq. (j3~4[) as a function of H. The dashed 
lines represent the roots v± in Eq. (35); the solid lines represent solutions of 
the cubic equation (|3~6l). 
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